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We consider a specific generalisation to spatial dimensions d greater than one of a formalism based 
upon conservation laws and the associated Ward identities, that exactly solves the one-dimensional 
Luttinger model, and expose its inadequacies for d > 1. We conclude that findings arrived at through 
application of this generalised formalism concerning systems of fermions with strong forward scattering, 
in particular that d = 1 would be the cross-over dimension from Luttinger- to Fermi-liquid state, are 
open to question. 
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I. INTRODUCTION 

The question with regard to the existence of non-Fermi hquid (NFL) metaUic states in spatial dimensions d greater 
than one, specifically in d = 2, has been subject of intensive investigations in the course of the past several years, 
following the observation that the normal-state properties of doped cuprate superconducting compounds markedly 
differ from those of conventional metals which are accurately described within the framework of the Fermi-liquid 
(FL) theory (Anderson 1988, 1989, 1990a,b,1997, Varma, Littlewood, Schmitt-Rink, Abrahams and Ruckenstein 1989, 
1990, Littlewood and Varma 1991). Since according to a celebrated theorem due to Luttinger (1961), validity of 
the many-body perturbation theory to all orders would imply that metallic states of fermions in spatial dimensions 
greater than one would be FLs, failure of the FL picture in the case of the doped cuprates has been considered 
to signal breakdown of the many-body perturbation theory as applied to these materials. Indeed, in the course of 
the recent years the question with regard to breakdown of the many-body perturbation theory has been a central 
aspect in debates concerning the above-mentioned failure of the FL phenomenology in describing the normal-state 
properties of the doped cuprate compounds (Anderson 1990b, 1993, Engelbrecht and Randeria 1990, Stamp 1993) (see 
also Engelbrecht and Randeria (1991), Anderson (1991)). Consequently, in this period a number of non-perturbative 
methods have been put forward and applied to models of interacting fermions with the aim of determining whether 
or not NFL states are viable in d > 1. Elsewhere (Farid 1999a, b) we have presented accounts of these developments 
and therefore refrain from repetition here. For our present considerations it is relevant to mention that the existence 
of NFL metallic states in realistic microscopic models of interacting fermions in d > 1 has remained elusive to this 
, . , date. 

In view of the above-indicated apparent discrepancy that prevails between theoretical predictions and experimental 
observations, we have undertaken to subject a number of aspects that are crucial to establishing feasibility or otherwise 
of NFL metallic states in d > 1 to an independent investigation. The present work is the third report of our 
investigations on this subject. Our first work (Farid 1999a) concerned a critical analysis of the above-indicated 
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^ Strictly speaking, this theorem has direct bearing on cases where fermions interact through a short-range two-body potential. 



1 



Luttinger (1961) theorem. We have demonstrated (Farid 1999a) that Luttinger's (1961) proof of this theorem involves 
an impUcit assumption which is strictly specific to FL metalhc states, implying that the Luttinger (1961) theorem 
far from excluding the possibility of the existence of NFL metallic states in d > 1 (for cases where the many-body 
perturbation theory is valid to all orders), in fact amounts to a statement of consistency of the FL states with the 
principles of the theory of interacting fermions as expressed in terms of the many-body perturbation theory. 

Our second work (Farid 1999b) concerned an in-depth analysis of the non-perturbative bosonisation scheme due to 
Haldane (1992) for metallic states in d > 1 as applied to a model of interacting fermions in d = 2, originally put forward 
by Houghton and Marston (1993) and subsequently extensively studied by Houghton, Kwon and Marston (1994), 
Houghton, Kwon, Marston and Shankar (1994) and Kwon, Houghton and Marston (1995). We have demonstrated 
failure of the Haldane (1992) bosonisation scheme in d > 1 and established that this scheme invariably predicts 
the metallic state of the model under consideration to be an unconventional FL, irrespective of the nature (and 
strength, provided it be non-vanishing) of the fermion-fermion interaction function, that is whether this be short- 
range, long-range Coulomb or super-long range. The latter result is in particular at variance with an earlier finding 
by Bares and Wen (1993), ^ indicating the metallic state of the system for super-long-range interactions to be a NFL. 
We have traced back the shortcoming of the Haldane (1992) bosonisation scheme in d > 1 to the neglect of some 
terms in the commutation relation concerning some current operators; in particular, in the case of metallic systems 
with a sufficiently small (or vanishing) value for the jump Z^^ in the momentum distribution function at the Fermi 
momentum, the neglected terms are more relevant than the sole term that is taken into account, namely the quantum 
anomaly. ^ We have also shown that under all conditions the Fermi energy of the interacting system within the 
Haldane (1992) bosonisation scheme coincides with that of the non-interacting system. This unambiguously implies 
that even for weakly-interacting metallic states, for which Zkj^ is close to unity, neglect of the indicated terms in the 
mentioned commutation relation for current operators is not justified. 

The exactly solvable one-dimensional Luttinger model (Luttinger 1963, Mattis and Lieb 1965) has been dealt with 
by three distinct methods. One is based on the technique of bosonisation, of which the Haldane (1992) scheme is 
a generalisation to higher dimensions. Q The second, due to Everts and Schulz (1974), is based on the solution of 
the equation of motion for the Green function which derives from an exact decoupling of the hierarchy of the Green 
functions. The third approach, due to Dzyaloshinskii and Larkin (DL) (1974), |^ exploits the conservation laws and 
the associated Ward identities. This scheme has been generalised to d > 1 and reported in a series of articles by Di 
Castro and Metzner (1991), Metzner and Di Castro (1993), Castellani, Di Castro and Metzner (1994), Castellani and 
Di Castro (1994) and extensively reviewed by Metzner, Castellani and Di Castro (MCDC) (1998). It is this scheme 
that we subject to a critical analysis in the present work. In our considerations that follow, we shall almost exclusively 
refer to the last review article without thereby implying that the cited subject matters would have their origin in this 
work. 

Contrary to case of d = 1, associated with linear energy dispersions for the non-interacting left- and right-movers 
(that is the case of the one-dimensional Luttinger (1963) model), in cases of d > 1 the generalised DL approach is 
not exact; MCDC (1998) assert, however, that for Fermi systems the generalised approach is asymptotically exact 
in the limit of strong forward scattering, that is in the limit of zero momentum transfer q for fermions in a narrow 
band of width A (A in the authors' notation) encompassing the Fermi surface. In this work we demonstrate this 
not to be the case. We further support this finding by referring to the similarity between the expressions for the 
single-particle Green function within the framework of the generalised DL (MCDC 1998) scheme and the Haldane 
(1992) bosonisation approach and the fact that, as we have indicated above, the latter is not adequate for d > 1 
(Farid 1999b). We also establish that the rigorous formalism for systems of interacting fermions, even in the limit 
of strong forward scattering, is not amenable to treatment by the Wilson-Fisher (1972) (Wilson 1973) technique of 



^ It is also at variance with the findings by Houghton, Kwon, Marston and Shanfcar (1994) and Kwon, Houghton and Marston 
(1995). We have however shown (Farid 1999b) that this discrepancy is due to use of certain invalid approximations by the 
latter authors. 

^ For completeness, we mention that by solely taking into account the quantum anomaly, the mentioned current operators 
form a Kac-Moody algebra (see, e.g., Goddard and Olive 1986). 

* The bosonisation method for calculation of correlation functions is due to Luther (1979) which derives from works by Luther 
and Peschel (1974), Mattis and Lieb (1965), Tomonaga (1950) and Bloch (1933, 1934). Parallel developments involving similar 
views as in the indicated works are due to Jordan (1935, 1936a, b), Kronig (1935a,b,c), Skyrme (1958, 1959, 1961a, b) and 
Coleman (1975). For reviews see Stone (1994) and Kopietz (1997). 

^ For a detailed exposition of this approach see (Bohr 1981). 
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dimensional continuation. 

The organisation of this work is as follows. In § II we present an outline of the generalised DL formalism (MCDC 
1998) which will facilitate our discussions in this paper. Anticipating the results that we present in § III, here we 
indicate that the conventional process of dimensional continuation is not applicable to a general formalism concerning 
systems of interacting fermions, not even in the limit where solely forward scattering processes are taken into account. 
In § III we analyse two crucial steps that are involved in the derivation of the generalised DL formalism and explicitly 
demonstrate their inadequacies in d > 1; in § III. A we deal with a relationship between the charge and the current 
part of a vertex function that has been held to be asymptotically correct in the limit of strong forward-scattering and 
show that its validity in d > 1 is strictly limited to the trivial case of non-interacting fermions; in § III.B we show 
the invalidity of a second relationship between two irreducible vertex functions, also held to be asymptotically correct 
in d > 1 in the limit of strong forward-scattering. In § IV we explicitly demonstrate how the singular behaviour of 
the momentum distribution function, in particular that pertaining to the ground-state of the system, at the Fermi 
momentum, renders the mathematical manipulations that are central to the derivation of the generalised DL method 
in d > 1 (MCDC 1998) inapplicable. Here we further show how similar manipulations, combined with a complete 
disregard of the errors associated with these, in d = 1 (specifically in dealing with the one-dimensional Luttinger model) 
give rise to results which coincide with the exact results, thus mislcadingly suggesting their general applicability even 
in d > 1. We finally present a summary and concluding remarks of this work in § V. We devote an Appendix which 
follows the main text to a brief exposition of the Wilson-Fisher (Wilson and Fisher 1972, Wilson 1973) technique of 
dimensional continuation, which concerns continuation of the number of the spatial dimensions from integer to real, 
as well as complex, values. From the details in this Appendix it becomes apparent how application of this technique 
is hampered in the context of the many-body theory of interacting fermions. Here we propose an alternative scheme 
of dimensional continuation which, although not unique, is not restrictive with regard to the functional dependence 
of the integrands on the variable of integration. 

II. OUTLINE OF THE GENERALISED DZYALOSHINSKII-LARKIN SCHEME 

Here we present a brief outline of the generalised DL scheme. For a detailed exposition of this scheme readers 
are referred to the article by MCDC (1998). In this Section we adopt the notational conventions of the latter 
work. For completeness, below p and q denote d+ 1-vectors (ipo,p) and (j(?o,q) respectively, with po (go) energy 
(throughout we assume h = 1) and p (q) momentum in d-dimensional momentum space (here d is integer; we shall 
be explicit when d is non-integer); {ipo, p)i, represents ipo when v = Q and — when ;/ = 1, 2, . . . , d; J^, /(fc), with 
f{k) EE i^(ifco;k), is the short for (27r)"i / dfco (271)"^^ / d'^fc F{iko;'k). For finite temperatures T, (27r)-i / dfco(...) 
symbolises a Matsubara sum, {2it)~^ J dko (...) ^ ^Sfco (■ " ^-'^^'^^ the summation runs over the discrete fermion 
Matsubara energies 27ri(fco -|- 1/2)T, with fcg = 0, ±1, . . .. Below, unless specifically indicated otherwise, we assume 
the summation convention A^B" = X)^!=o ^vB" . 

For a system whose Hamiltonian H is of the form 

H = Ho + Hj-.^J^^l^Ka^K'^ + ^ 3kk'(q) «Lq/2.<T«k+q/2,<TaL+q/2,'T'"k'-q/2,^', (1) 

k k,k',q (7(7' 

with dj^ ^, dk.<j the canonical creation and annihilation operators for fermions with spin a, respectively, the energy 
dispersion of the non-interacting fermions, ^ V volume of the system and ^^k' (l) fcrmion-fcrmion interaction 
function, MCDC (1998) attempt to transform the expression for the self-energy Sct(p), 

into a closed form. Here Dpp,{p ~ p') is the screened interaction function (see Eq. ([l^ ) below), Ga{k) = Ga{iko;h) 
the (thermal) single-particle Green function. 



^ We assume = ||k||^/2 (in units where nif, the fermion mass, is equal to unity). For some relevant remarks see Footnote 
id further on. 
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A°'^(p;9):=(n<T'('?)ap-«/2,<Tap+g/2,<T)ir'' (3) 
the irreducible density (or charge) vertex, in which |^ 

"<T(q):=Xl^k-q/2,<T«k+q/2,'^- (4) 

k 

In Eq. "tr" indicates that the external Green functions are not taken into account (i.e. they are truncated), 
which amounts to dividing the expression associated with the (thermal) average (. . .}'"', without subscript "tr", by 
Ga{p+q/'2)Ga-{p—q/'2); the superscript "irr" indicates that in the diagrammatic expansion of the associated correlation 
function only one-interaction-irreduciWe diagrams are taken into account, i.e. those which do not fall into two disjoint 
parts upon removing an interaction line in them. Through reliance upon the generalised closed-loop theorem (CLT) 
(Kopietz, Hermisson and Schonhammer 1995), |^ MCDC (1998) relate A'^,^{p;q) with the irreducible vertex part 
associated with the noji-mteractiiig density-current jgiq), that is 

A^(P; 9):=(jo (9) "■p-q/2,>7a-l+q/2jtr' (5) 

as follows: 

Al,Jp;q)^S,,,,Al{p;q), A,, ,(p; g) = ,5,,,, (p; g) . (6) 

Here the d-dimensional vector A^i „{p] q) corresponds to the u — 1,2^ . . . ,d components of AJ^, ^(p; q). The expressions 
in Eq. (^ can be written as 

AO ^ A r (igo,q)^A^(p; q) 
«9o - q • 

The density-current vertex function A[^ (p; q) is related to the vertex part 

KiPi q)--Hf {q)ap-q/2,aal^q/2,a)^^ (8) 

and the irreducible vertex part corresponding to the interacting density-current j'^{q), that is 

A/(p; g):=(?(q) a,_q/2,aal^q/2,X , (9) 

through the Dyson equation 

r:;(p;<z) = A/(p;g) + i^ J] r'"^' {q) g^, ^Z, (q) A^" ^" {p; q) , (10) 

a' k" a" 



where 

r'"'{q):=- j ^ (rl^K+^AA-.A.) ^ / '^Tl{k-,q)G,{k + q/2)G^{k-q/2), (11) 

(«Lg/2,^'afc+9/2,-'ap-9/2.^a|;+5/2,a)i7- (12) 

Owing to the exact Ward identity 

(z9o,q).r^(fc; q) = G-\k + q/2) - G-\k - <7/2), (13) 



Note that y^^nj(g) = jo{q) = fil): where jo{q) and j'^{q) stand for the density (or charge) parts of jo{q) and j'^{q), 
respectively (see Eqs. ( p^ and ( ^ ) below). 

* The CLT (or, loop-cancellation theorem) for the one-dimensional Luttinger (1963) model is due to Dzyaloshinskii and Larkin 
(1974). For an extensive exposition of this theorem see Bohr (1981). 
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from Eq. (11) it readily follows that 

{tqo, ci),r-^%q) = n,(k - q/2) - n,(k + q/2), (14) 

where 



G^{k±q/2) 



(15) 



stands for the momentum distribution function corresponding to the spin-cr fermions. Through relying upon (see 
however §§ III and IV) n^(k ± q/2) = n^(k) ± iq • Vkn^(k) + Odlqf), for smaU ||q||, MCDC (1998) deduce from 
Eq. (|l|) (c.f. Eq. ^ below) 

{^qo,q).r■'^^{q) = -q • Vkn.(k) + 0(||qf ). (16) 
The screened interaction D^^,{q) in Eq. (||) is that within the random-phase approximation (RPA), 



D^i>{q) = g££;(q) f .9££;;(q) G{k" q^ak" + q/2)Di:,i,{q). 



(17) 



In practice, the interacting 'bubble' G{k" — q/2)G{k" + q/2) is replaced by the non-interacting one (MCDC 1998), 
i.e. Go{k" — q/2)Go{k" + q/2); this on account of the generalised CLT, under the assumption that k" is confined to 
a small neighbourhood of kj-. For a comment on the generalised CLT in d > 1 see (Farid 1999b). 
From the equation of continuity [n(q),iJ]_ = q • j(q), where 

n(q):=.^n.(q) (18) 

cr 

(see Eq. (^ above), one can obtain the conserved current j(q) for the system governed by the Hamiltonian H in 
Eq. (p. The conserved current jo(q) corresponding to Hq in Eq. (|l|) is obtained through application of the continuity 
equation [n(q), 7?o]- =q-jo(q),Q 

jo(q) = 51 ak-q/2,<T«k+q/2,<T, (19) 

where v|^:=Vk£k ^ ^^(^ commutator [n(q),iJ/]_ gives rise to 'non-diagonal' quartic contributions (involving four 
different momenta) to the current operator j/(q):=j(q) — jo(q), which are irrelevant in the low-energy limit (from the 
renormalisation-group perspective); taking the 'diagonal' contributions, corresponding to the momentum transfer q' 
equal to either q or — q, the following (approximate) expression is obtained (MCDC 1998) 



[n{q),Hi]_ = ^J2J2 [fi.(k-q/2)-n,(k + q/2)] 4,^^/2 ,,ak'-hq/2,.' 

k fj k' fj' 

+ ^ E E «Lq/2,>k+q/2- [fi-' (k' - q/2) - n,' (k' + q/2)] . (20) 



2V 

k (T k' cr 



MCDC (1998) cast the right-hand side (RHS) of this expression into the form q • j/(q), with 

j/(q) = E '^k^(q) «k-q/2,cr«k+q/2,<T, (21) 



where 



^ Explicit calculation yields [n(q),i5'o]_ = J2ka{^i+i/2 " Sk-q/2) aLq/2,<T^k+q/2,<T- With el = ||k|p/2 (see Footnote 
above and Footnote |l^ below), one has £k+q/2 ^ Sk-q/2 = q • k = q ■ v[^. 
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-'L(q):=--^E ff£k''(q)Vk'n..(k'). (22) 

k' CT' 



In obtaining the result in Eq. (|2l|), MCDC (1998) employ n<,(k - q/2) - n^(k + q/2) = q • [Vkn<T(k)] + 0(||qf ) 
(see §§ III.B and IV however; also compare with Eq. ( |6^ ) below), and the same for ncr'(k' — q/2) — no-'(k' + q/2). 
On account of the fact that for the ground and low- lying excited states of the system [Vkncr(k)] is most significant 
for k = ki? (the same concerning [Vk'n(j'(k')] for k' = k^), MCDC (1998) make use of the following relationship 
concerning the radial part of the k' integration (that is that in the direction normal to the Fermi surface) in the 
expression for v^^(q) in Eq. (^2|), P"| 

J dfc; fc;'-'5££''(q) Vk'n.'(k') « ||k^r-\g££;^(q) J dfc; Vk'n.,(k') 

= -\\k'p\f-'g^^,Jci)n^,^, (23) 

where nk'^ denotes the outward unit vector normal to the Fermi surface at k'^. Here d is not necessarily integer 
(see Appendix). In arriving at the last expression on the RHS of Eq. (^), MCDC (1998) have made use of the fact 
that for the ground and low-lying excited states of the system, n(j'(k') yields unity for k' deep inside and zero for 
k' far outside the Fermi sea (that is for k' = k^ and k' = k^, respectively). MCDC (1998) observe that the RHS of 
Eq. ( p3| ) (which is a c-number) would remain intact if no-'(k') on the left-hand side (LHS) were the expectation value 
(n„'(k')) = n„'(k'). Q Thus MCDC (1998) employ (c.f. Eq. (||) above) 

vL(q) - vL(q):- - E 3k£''(q) Vk'n,,(k'), (24) 

k' a' 

and consequently (c.f. Eq. ( ^ above) 

j7(q) ^j/(q):=E vL(q)aLq/2,^«k+q/2,^- (25) 

ko- 

It is this j/(q) that MCDC (1998) employ in the construction of the generalised DL scheme. In § III.B we extensively 
discuss the shortcomings of the foundation of this expression for j/(q). Here we only briefly mention that, contrary to 
the statement by MCDC (1998), the relationship between the second and the third expression in Eq. ( p3| ) is defective: 
momenta being confined to a narrow band of width A circumscribing the Fermi surface, we have ||ko — ki|| = A. 
With < A <C ||ki?||, unless the ground and low- lying excited states of the system be uncorrelated (or at most 
weakly correlated), it is not permitted to replace no-'(ki) by 1 and no-'(ko) by 0. To clarify this point, let us take the 
expectation value of the expressions in Eq. ( p3| ) with respect to the A^-particle ground state of H. With < A ^ Ilk'p-H, 
it is seen that —1 on the RHS of Eq. ( |2^ ) has to be replaced by r\a-'(ko) — na-'Cki) « — Zk'^o-': where Z]f^i^„i stands for 
the amount of discontinuity in <n„i{]<i') at k' = k'^ which is less than unity for interacting systems and is equal to zero 
for such NFL states as the marginal (Varma, et al., 1989, 1990, Littlcwood and Varma 1991) and Luttinger liquids 
(Haldane 1980, 1981, Anderson 1997). This undue replacement of ^k'^^o-' by unity, is exactly the step undertaken in 
the construction of the Haldane (1992) bosonisation scheme in d > 1; the consequences of this replacement would not 
be as severe were it not that in cases where Zk^ a' is small, or vanishing, the other contributions, which are neglected 
within the bosonisation scheme, were not so significantly more dominant than that of the quantum anomaly that in 
an enhanced way (through replacing Zk^ a' by unity) is taken into account (Farid 1999b). We observe that, unless the 
system under consideration be non-interacting, the constraint on A, namely A <^ ||ki?||, implies that in contradiction 
to the suggestion invoked by Eq. (p3|), / dfc'r fc'^^^ .g^^/ (q)Vk'no-'(k') is not even approximately a c-number. Be it 
as it may, as we shall demonstrate in § III.B, the expression in Eq. ( [2l| ) with v^^(q) as defined in Eq. ( p2[ ) is not the 
appropriate expression for j/(q) and, moreover, the appropriate expression (see Eqs. ( |5l| ) and ( p9| ) below) even for 
large A cannot be replaced by a c-number. 

Through applying {iqo,q)„ to both sides of Eq. ^ and making use of Eqs. (|l|), (|l|) and (|2|), MCDC (1998) 
obtain the following crucial relationship 



^'^ In the thermodynamic limit one has V ^ X/k' (•••)"* (27i") / d'^k' (...). For d-dimensional integrations see Appendix. 
MCDC (1998) state namely "The k'j. integral over the operator expression thus yields a c-number, and the same c-number 
is obtained when replacing h^i(k') by its expectation value." 
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{iqo,q)^r'^{p;q) = {iqo, q)„A'^{p; q). 



(26) 



In arriving at this result, use has been made of the relations (27r) ^ / dfcg aj./_^^2 cr'°'fe'+9/2,o"' = '^k'-q/2 cr''^k'+q/2,CT' 
and 

- KiP', Q) + {ii{q)ap-q/2,aal^g/2jtJ = (jo (9)ap-g/2,aap+g/2,o->i7 ^<t(p; q)- (27) 
From the exact Ward identity in Eq. ( p^ ) one finally obtains 

(zgo, q).A^(p; q) = G;i(p + g/2) - G'^p - q/2) (28) 
which in combination with Eq. ^ yields (MCDC 1998) 

Eq. (||) in conjunction with Eqs. (p7|), (|2^) and the Dyson equation leads to a closed set of equations for the 
single-particle Green function Ga{p)- One has 



Ga{p) = Go;.(p)|l - / ^PP^(P P) G„{p')+X,{p)G„{p)\, (30) 

1^ Jp' i{p - Po) - V[p+p,]/2 • (P - P') J 



where Go-^ip) denotes the non- interacting counterpart of Ga-{p) and 



Jp' z(p-po)-v{'p^p,j/2-(p-p') 



In practice (MCDC 1998) D^^,{p-p') (see Eq. 1^)) is approximated by D^^{p - p')^:D^'' (p - p'). 

For the case of an isotropic system, through replacing ||v|^|| by Vp and employing for small q:=p — p' (c.f. the 
denominators of the expressions on the RHSs of Eqs. (pG) and (pl|)) 



vfp+p']/2 ■ (P - P') = VpiPr - P'r) + 0(l|qf ), (32) 

where Pr:=||p|l - ilkF||, Pr-=llP'll ^ II^fII, and for small q = {iqo, q), 

D;^q)^Diqr/qo,qt/qo), (33) 

where qr-=<l ■ p/||p|| and gt:=(||q||^ — q'^Y^'^ > 0, the radial and tangential components of q, respectively, MCDC 
(1998) render the expression for Go-(p) in Eq. (|30| ) suitable for treatment by the conventional method of dimensional 
continuation (see Appendix — we note in passing that here qt is the component of q inside the tangent space specified 
in Appendix). In this way MCDC (1998) establish that the cross-over dimension from Luttinger- to FL state is d = 1; 
that is, the metallic state of the system under consideration is a FL (though not necessarily a conventional one) in 
d > 1. This finding is in agreement with the e-expansion result by Ueda and Rice (1984) around d — 1; according 
to these authors, in the d = 1 + e fermion system the only weak-coupling fixed point is the FL one. As we shall 
demonstrate in § III. A, the expression in Eq. (E9h whose denominator on the RHS is the origin of the denominators on 



the RHSs of the expressions in Eqs. ( |30| ) and ( pl[ ) (which is further simplified in Eq. (32)), is not valid except in the 
strictfy non-interacting case. In § III. A we further indicate that for NFLs, I]CT(p±q/2,cj±cjo/2) does not allow for an 
effective disentangling of q from p (as well as luq from to) and consequently for these systems the appropriate equation 
for Gcr{p) does not accommodate terms involving contributions whose form would resemble the denominator of the 
integrands on the RHSs of Eqs. ( |30| ) and (|T]) which thus would be amenable to such reduction as in Eq. (^). Note 
that the expression in Eq. ( ^2| ) is specific in that it does not involve a component of p' in the tangent subspace S± 
(see Appendix). Thus we conclude that the general formalism for systems of interacting fermions even in the limit of 
strong forward scattering does not lend itself for treatment by the conventional method of dimensional continuation. 

In the following Section we subject the generalised DL formalism to a critical analysis and demonstrate its failure 
in d > 1. We also indicate the underlying reasons for the exactness of this formalism in d = 1 when applied to the 
one-dimensional Luttinger model (1963). 
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III. CRITICAL ANALYSES 



In this Section we present our analyses of some crucial elements that are involved in the derivation of the generalised 
DL method described in § II. 

A. On the relationship between A°(p;g) and Xa{p;q) 

Here we restrict our considerations to the case of zero temperature (T = 0) and deal with the analytic continuations 
of functions of imaginary energy along the real energy axis. We denote the analytic continuations of such functions 
as Ga{p) and T,c,{p) by Gc{p,uj) and T,c{p,uj) respectively. 

From the Dyson equation for the Green function it follows that fn 



G-\p + q/2, u; + LOo/2) - G'^p - q/2, lu - loo/2) = t^o - q • 





per 



- [S,(p + q/2, LO + LUo/2) - S,(p - q/2, lo - uJo/2)] . (34) 

Assuming continuous differentiability of S(j(p', uj') with respect to p' and w' in neighbourhoods of p' = p and ui' — ui 
respectively, for sufficiently small ||q|| and |wo| we have 

G-\p + q/2, u + iOo/2) - G-\p - q/2, u; - loo/2) = ujo [l - 9I],(p, u;)/dLo] 

-q- [vO, + VpI].(p,u;)] +O(||qr,|t^or'), (35) 

where 7,7' > 1. From Eqs. (^) and (|3|), up to corrections of the form C'(||q||''', jcjol'*' ), we deduce 



A°(p,c.) = l-^^^^, A.(p,c.)-vO, + VpE.(p,c.). (36) 
Following the expressions in Eq. m) relating A" to A^-, the results in Eq. (pq) imply 



V^^Ap,co) = - ^^^^:-' v°„ (37) 



dT,^{p,Lu) 
doj 

which must hold exactly for p = p^ and uj — uip the results in Eqs. (^ and (^) are to be asymptotically valid 
for ||q|| and \ujo\ — > 0. We note that a T,cr{p,uj) whose associated T,cr{pF,t^) and I]ct(p,wf) are continuously 
differentiable functions of uj and p in neighbourhoods of uj = lof and p = P_f, respectively, by definition corresponds 
to a FL (Farid 1999a). In other words, if valid, Eq. ( |37| ) must hold for FLs when lu ~ ujp asp^p^ and p = P_f as 

LU — > LUp. 

Let us now consider the following asymptotic expressions which apply to conventional FLs in rf = 2 for p close to 
the Fermi surface and lu —^ luf (Hodges, Smith and Wilkins 1971, Bloom 1975) (in what follows we assume the 
system under consideration to be in the paramagnetic phase and thus suppress the spin indices): 

ReJ:ip,Lu)^j:ip,LUF) + f3p{uu-LUF), (38) 
ImS(p, Lu) ~ ttp sgn(w — LUp) (^^ — wf)^ In \uj — ujf\, (39) 

where Pp and Ofp are real valued and for interacting FLs /3p < (strictly negative; see text following Eq . (|4^ ) below) 
and Qfp > 0; for conventional FLs in c? = 2, E(p, w^) ^ and ap > 0. Substitution of the RHSs of Eqs. (|38|) and <^ 
in Eq. (p7[) results in 



12 In expressions akin to that in Eq. (|^ (such as that presented in Eq. Q below), MCDC (1998) add OdlqlP) to the RHSs 
in order to account for the fact that in general the non-interacting energy dispersions can be more complicated than linear or 
quadratic functions of momentum. We do not follow this practice for the reason that the one-particle kinetic-energy operator 
Ho and the velocity operator v are related to one another through the momentum operator p; we have Ho-~p^ /2 and v:=p 
(we use mf — 1 throughout), so that in the momentum representation must hold (assuming a uniform state) £p = ||p||^/2 and 
Vp = p for all p (from this perspective ~ Sp-q/2 = 1 ' ^pi whence absence of OdlqlP) in Eq. (^2|)). A fully consistent 

generalisation of the formalism considered in the present work, that is one capable of dealing with general non-interacting 
energy dispersions, must therefore involve some reformulation of the problem at hand. 
See also Fujimoto (1990) and Fukuyama, Narikio and Hasegawa (1991). 
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VpS(p, ujf) + [Vp/3p] {uj - cup) ~ -/3p vO , (40) 
[VpQfp] {u;~u;f) \n\uj - ujf\ (2 In |cj - cjf| + l) ap v° . (41) 

We now consider the consequences of these resuhs. As will become evident, Eqs. (|40| ) and (^l|) hold only for strictly 
non-interacting FLs (in this case they hold identically), thus implying incorrectness of Eq. (^7|) and by extension that 
of the expressions in Eq. (^) for interacting systems. For lu = ujf, Eq. (|40| ) leads to 'VpT,{p, ujf) = ^Pp'^p (for p in 
the vicinity of the Fermi momentum), which, according to Eq. (^ ) ior u; luf, implies Vp/3p = /3p = (see 



further on). This result in turn necessitates Wp'S{p,luf) = 0; further, for uj — > ujf the LHS of Eq. (41) vanishes while, 
unless Qfp = 0, the RHS diverges. It follows that unless the system of fermions (in d = 2) be non-interacting, Eq. ( |37| ) 
and, in consequence, the relations in Eq. (||) are incorrect (see text following Eq. ( ^ ) below). For non- interacting 
systems we have namely (see Footnote [ij above) 

Go\p + q/2, u; + u;o/2) - G^^p - q/2, u; - luo/2) = ujo ~ q- (42) 

from which it is readily deduced that, with (IjVp)" the non- interacting counterpart of the irreducible vertex part 
A[^(p;g), the relations in Eq. (|^) are identically valid in the strictly non- interacting case. This completes the proof 
of our statement with regard to the condition for validity of the expressions in Eq. (^. We emphasise that the case 
of d = 1 is very special when Vp is a constant, say v*^ (which points to a definite direction), as is the case for the 
one-dimensional Luttinger (1963) model. Here the relations in Eq. (^) hold identically, for whatever A° (p; q) may be, 
from the definition for jg (q), namely (see Eqs. (p^), (^) and ( p^ ) above) 

Jo(q):-(^(q),jo(q))'' ^E(l,v°)^4_q/,^ak+q/2., (43) 

ko- 

it follows that A^{p]q) is v*^ times A°(p;g) (see Eq. (^)). 

We conclude that the relations in Eq. (^ are not asymptotically correct in the forward-scattering limit when d > 1, 
unless the system under consideration be non-interacting, in which case they are identically valid in any spatial 
dimension. 

Above we have asserted that Vp/?p = ^==^ /3p = 0. We now demonstrate the validity of this equivalence. For 
FLs, it can be shown that (Farid 1999a) 



1 f°° du' 

f3p = - —2 [lmS(p, u;' + u;f) ~ ImS(p, -lj' + u;f)] ■ (44) 
^ Jo ^ 

We note in passing that since for interacting FLs the integrand of this expression is not identically vanishing and, 
moreover, owing to the assumed stability of the system lmY,{p,uj' + lof) and — Iml](p, —uj' + lof) are never positive 
for uj' > 0, it follows that /3p < (strictly negative). The expression in Eq. (^) can be transformed into |^ 

1 dz 

f3p^-lm [E(p, z) - E(p, 2(1 - z)] , (45) 

where /j, S {^n , fJ-N+i) with ^j\j:^Ej\jfi — En^i^q = ujf and ^n+i'-—En+i.o — Ej\jfi. Here Ema stands for the total 
energy of the interacting M-particle ground state and N denotes the actual number of the particles in the system. The 
2;-integration in Eq. ( p5| ) is carried out along an arbitrary contour in the upper half of the complex z plane, connecting 
2 = yU -|- 0"*" with the point of infinity. The arbitrariness (under the mentioned constraints) of the contour of the z- 
integration has a far-reaching consequence: assuming differentiability of /3p with respect to p (in some neighbourhood 
of Pf), it follows from Vp/3p = that Vp [E(p, z) — S(p, 2/i — z)] =0 for all z in the upper half of the complex z 
plane; by the Riemann-Schwarz reflection property with respect to the real axis, this identity holds also for all points 
in the lower half of the complex z plane. |^ This implies that E(p, z) = S(p, 2/i — z) over the entire complex plane. 



For FLs dz (z — fi) ^ Im[E(p, z) — E(p, 2/i — z)] exists, however J^^°° dz {z — fi) [S(p, z) — E(p, 2/i — z)] does not. 

Addition of the infinitesimal 0"*" to /t in the lower boundary of the latter integral (see Eq. (^^) formally renders this integral 
existent without affecting the value of /3p. 

For Im(z) 7^ 0, E(p, z) satisfies the Riemann-Schwarz refiection property E(p, 2*) = S*(p,2;) (see e.g., DuBois 1959, 
Luttinger 1961, Farid 1999a). 



9 



In view of the expressions in Eqs. ( [3^ ) and (|39|), this result is seen to apply only if E(p, 2:) is independent of ^, a 
property which is specific to noiz-interacting systems. A corollary to this result is that /3p = (that is zero rather 
than any other constant that Vp/3p = would suggest in the first instance; see the text following Eq. (El|) above ). p| 
Our above analysis has explicitly concerned FLs (through our explicit use of the relations in Eqs. ( |38| ) and (|3S|), 
but more importantly Eq. (^4|)), however, as can be easily verified, for NFLs the RHS of Eq. (^^ does not allow for 
any effective separation of q and ujq from p and to respectively, so that in > 1 the relations in Eq. lack even the 
most fundamental basis for validity as regards NFLs. 



B. On the relationship between A'^{p;q) and r^(p;g) 



Evidently for obtaining the j/(q) as presented in Eq. (|2l| ) in terms of v^o-(*l) ^'i- dH); the functions in the square 
brackets on the RHS of Eq. (20) have been expanded to linear order in q. Such an expansion is in general invalid for 
reasons that we discuss in some detail in § IV. A second step involved in the derivation of j/(q) in Eq. ( p5| ) is that 
of replacing VknCT(k) by the expectation value Vk;(ncr(k)). The justification that has been given by MCDC (1998) 
for this substitution is untenable (in particular in cZ > 1) (see also text following Eq. (|2^) above). For obtaining a 
generally valid expression for j7(q), the RHS of Eq. (|20| ) must be first written in terms of four separate contributions, 
each involving only one of ncr(k — q/2), no-(k + q/2), no.'(k' — q/2) and na'{^' + q/2). Subsequently, in each of these 
contributions, by appropriate transformations of the summation variables k and k', the vector q in the arguments 
of these functions must be transferred into that of the interaction function 5kk'("F^)' which we assume to be a 
continuously differcntiablc function of its indices k and k' (c.f. statement preceding Eq. (7) in Castellani, Di Castro 
and Metzner 1994). In this way one obtains (for small ||q||) 



j/(q) = ^J2Y1 {['^k.9££''(-q)] n.(k)aL_q/2,.'ak'+q/2,^' 

k o" k' o"' 

+ [Vk'3k£''(q)] «Lq/2,.«k+q/2,. n.'(k')}. (46) 

Making use of the commutation relation 

^k-q/2,(T"k+q/2,o-, no-'(k') ^ = Sa^a' {^k+q/2,k' " '^k-q/2,k' } ak_q/2,o- ak+q/2,o-, (47) 

which is obtained through a repeated application of the canonical anti-commutation relations for aj^. ^ and ak,o-, the 
expression in Eq. ( ^6[ ) can be written as 

j^q) = ^k^(q) "Lq/2,<T«k+q/2.<T 

kcr 

+ ^E{['^k' 5£^'(q)lk'=k+q/2] - ['^k' 5^^'(q)lk' = k-q/2]}aLq/2..«k+q/2,., (48) 
k a 

where (c.f. Eq. ( [2^ ) above) 

vL(q):-^ E ['^k'ff£^:(q)] fi.'(k'), (49) 



k' a' 

in which 



5££'(q):=^(ff££'(q) +ff£''£(-q)). (50) 



Let Re(z)=:(jj > /i and Im{z) — :ri, rj I 0. For such a z we have ImE(p, z) < while ImE(p, 2/i — z) > (validity of both of 
these relations is necessary for the stability of the system). Satisfaction of E(p, z) = E(p, 2/i — z), for all z, in this particular 
case would imply ImE(p,tj + irj) = for all u > /j.; it similarly would imply ImE(p,ct; — irj) = for all u) < fi. By the 
Kramers-Kronig relation it follows that E(p, 2) is independent of z. Note that in obtaining this result we have not relied upon 
the specific form of the expressions in Eqs. (k8) and (|39). 
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By assuming Vk'.gkk'(Q) to be a continuously-difFerentiable function of k', the last term on the RHS of Eq. ( p^ is 
seen to be proportional to q, so that in the limit of q — > it can be neglected in comparison with the first term. For 
small ||q|| we can thus write 

= Yj ^k^(q) "I-q/2,^^k+q/2,<T, (51) 
kcr 

where v^^(q) is defined in Eq. (^). Note in passing, that for a constant interaction (i.e. a k, k'-independent g^^, (q)), 
VkCT(q) = 0, so that Eq. (1^) correctly yields j/(q) = 6, or j(q) = jo(q). 

The correlated velocity operator in Eq. (^ ) is distinct from that in Eq. (^2|) . To bring out the fundamental difference 
between the two, let us follow MCDC (1998) and replace the operator fio-' (kMon the RHS of Eq. (|49| ) by its expectation 
value n(j'(k') and compare the resulting vecto r v£ „(q) with that in Eq. (pij). Since by assumption Vk'gkk'('l) ^ 
smooth function of k' (see text preceding Eq. (|4q ) above), we can write (recall that — k^|| = ^ ^ ll^^ll) 

J dK k'/^' [Vk'5££''(q)]n.'(k') « [ Vk'5££:(q)|^,^^, ] J dK k'/-' n,,(k'). (52) 

For the non-interacting case, where no.'(k') = 0(k'^ — k') (recall that distances are measured along ^k^), the in- 
tegral on the RHS of Eq. (||) is equal to (|ik^||'^ - ||k^||'^)/d, so that the RHS of Eq. (||) becomes equal to 
[^k'5££''(q)lk'=ki,] (Ilk^.H'^ - ||k^||'')/d, to be compared with the RHS of Eq. (H), which is equally applicable to 
the non- interacting ground states, namely — 5kk' (l) ll'^i^ ll'' '^k^ ■ The two results are manifestly different. This 
demonstrates that Eq. (|2j), and therefore Eq. (^5|), is fundamentally incorrect. This aspect is additional to the 
shortcoming of the result in Eq. (|2^), which as we have discussed in the text following Eq. (^5|), shows up when the 
ground and low-lying excited states of the system are (strongly) correlated. Thus we have unequivocally established 
that not only the change in Eq. (p5|), corresponding to replacing v^cr(ct) as defined in Eq. ( |2^ ) by v^^(q) as defined in 
Eq. (|2j), lacks justification (unless the system under consideration be non- interacting) , but also that the expression 
for v^(q) as presented in Eq. ( p^ ) (to be compared with the v^c,(q) in Eq. (^)) is invalid, this on account of the 
fact that it relies on the Taylor expansion of an operator (that is no-(k±q/2)) whose, in particular, expectation values 
with respect to ground states are in general non-differentiable (see § IV). 

A major consequence of use of Eq. (|T]), or Eq. ([2l|), rather than one in which the operator ^ (q) is replaced 
by the c-number v^^(q) (effected through substituting nCT'(k') for ncr'(k') in Eq. (^), or Eq. (p2[)), is violation of 
Eq. (|2^) which is the crucial link between the exact Ward identity in Eq. (|l^) and the expression in Eq. (|2^). To make 
this statement explicit, we first apply (iqo, q)iy to both sides of Eq. ( p^ ) and subsequently substitute the expression 
on the RHS of Eq. (|l|) for {iqo,q)^,r'^' (q). Note that in doing so, we depart from the approach by MCDC (1998) 
who employ the expression on the RHS of Eq. (p^). Next, we write the expression for {iqQ,q)^J'^'^ °' (q) in the 
contribution corresponding to {iqQ,q)^T'^{p;q) as two separate terms and in each of these by transformations of the 
dummy variables of summation, remove q/2 from arguments of n(j'(k' — q/2) and n^'Ck' + q/2). This procedure is 
exactly the same as that which we employed for deducing j/(q) from the equation of continuity associated with Hj 
(see Eq. (Eq) above and text following it). Thus we obtain 



E im,cD..r''''^\q)9^:i:{q)A^"-"{p;q) 

E {-^E"-'(l^')['^k'5£''£"(q)]}Ar-"(p;g). (53) 



V 

k' a' \l" U' 



k" a" k' <7' 



Neglecting the possible difference between 5kk''(q) ^^^^ 5k k' (q) i^'^^ ^'i- ( |50| ) above), it is seen that the term enclosed 
by curly braces in Eq. (^3|) is equal to the expectation value of ^ (q) as presented in Eq. ( ]49| ) . As we have extensively 
discussed above, there exists no justification, whatever, for v^^(q) to be a good representation of v^o.(q). If, as in 
the treatment by MCDC (1998), 'v^^{q) could be replaced by its expectation value, then making use of the defining 
expression for A^" (p; q) in Eq. (|l|) and the expression for j/(q) in Eq. ( |5l] ) wherein "v^^i^q) is replaced by v^^(q). 



the RHS of Eq. (53) would reduce to q- Aa-{p; q). Adding this result to (iqo, q)yA^'^(p; q) (see Eq. (^0|), to whose both 
sides {iqo, q)i/ is applied), one would arrive at Eq. (p6[). We emphasise that in arriving at Eq. ( |2^ ) it is the substitution 
of '^k cr (q) ^k cr (q) ^^^^ plays the vital role and not whether ^ (q) is the expectation value of ^ (q) as defined 
in Eq. ( p^ or Eq. (^). However, it should be realised that contrary to the case of the v^jj.(q) in Eq. ( |2^ ) whose 
replacement by its expectation value can be justified for non-interacting systems (and by extension, weakly interacting 
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systems) (see § II, text following Eq. (p5|)), replacement of the Vko-Cl) (49), which we have demonstrated to 

be the correct velocity operator, by its expectation value is not justified, not even for cases where the ground and 
low- lying excited states of the system are non-interacting (or weakly interacting). Of course, in the non- interacting 
case ffkk'(q) identically vanishing so that the problem associated with the inadmissibility of replacing v^^(q) by 
v^^(q) does not arise in this hmit (see Eq. ( ^ ) above). 

Above we have demonstrated that in d > 1, the trivial relationship between the vertex parts A[^(p; q) and T'^{p] q) 
as given in Eq. (^) is invalid and consequently no direct link exists between the exact Ward identity in Eq. ( p^ ) and 
the expression in Eq. (p8|). Further, our considerations in § III. A made evident that also Eq. (^ has no rightful place 
in the theory of interacting fermions in d > 1. 



IV. MATHEMATICAL TRANSGRESSIONS AND "BLESSINGS" OF INCORRECT ERROR 

ESTIMATION 

In this Section wc illustrate how in d = 1, specifically in dealing with the one-dimensional Luttinger (1963) model 
at r = 0, a combination of transgressions against mathematical rules and uncritical neglect of truncation errors can 
give rise to outcomes which coincide with the exact results, thus creating the false impression that this approach 
were justified. The particular example that we consider in this Section makes explicit that specifically in d > 1 the 
same procedure yields manifestly incorrect results. Our illustrative example involves the ground-state momentum 
distribution function n(k) which is singular at k = ki? both in the case of FLs and NFLs. 

We first explicitly deal with the one-dimensional Luttinger model for spin-less fermions (Luttinger 1963, Mattis and 
Lieb 1965), characterised by the anomalous dimension a, a G (0, 1). We note in passing that for the one-dimensional 
Hubbard model < a < 1/8, with a = 1/8 corresponding to the limit of infinite on-site repulsion (Voit 1993). By 
choosing the positive direction to be from left to right, below we consider the one-dimensional vectors k and q as 
scalars and kp = Ik^l- Here the ground-state momentum-distribution function n,.(k), with r = ± the branch index 
(r = + indicates the branch of the right movers and r — — that of the left movers), is continuous but not differentiable 
at k = rkp. For the density-density correlation function Hrril) pertaining to the interacting system we have (DL 
1974; see, e.g., Metzner and Di Castro 1993, Voit 1994) 

1 a (a) /.A/2+rfeF-q/2 

n..(g) ^ — . ^''^ 0^ , 5.(q):= / dk [n,(k-q/2)-n,(k + q/2)], (54) 

27r iqo - Vj! • q J -X/2+rkF+q/2 

where vjl stands for the constant velocity of the non-interacting fermions; it points to the positive direction for r = + 
and the negative direction for r — ^. Further, A is a cut-off which is finite but can be arbitrarily large; for our 
following considerations it is required that A be sufficiently large so that nr(k) can be replaced by zero (unity) for 
r = + (r = — ) when k « A/2. The freedom to chose A arbitrarily large is a specific aspect of the one-dimensional 
Luttinger model. We note that although introduction of such a cut-off is necessary for a correct treatment of the 
Luttinger model (Mattis and Lieb 1965), imposition of a 'sharp' cut-off on nr(k), implying n±(=FA/2 =F 0+) = and 
n±(^A/2 ± 0+) = 1, gives rise to an incorrect asymptotic result for the long-distance behaviour of in particular the 
single-particle Green functions Gr pertaining to the model under consideration; |^ for this long-distance asymptotic 
behaviour to be described correctly, it is necessary that n+(— k) and n_(k) are reduced to zero for k ^ A/2, which 
can be achieved through imposing a 'soft' cut-off on these functions which becomes operative for |k| larger than A/2 
(Farid 1999b). 

Through decomposing the integral in Eq. ( p^ into two integrals involving separately n,.(k — q/2) and n,.(k -|- q/2) 
and subsequently transforming the integration variable k, with k ^ k + q/2 in the first and k ^ k — q/2 in the second 

integral, cancellation of a common contribution to the resulting integrals gives rise to 5r(q) = I-\/2+rkF'^'^ nr(k) — 
I\/2+rkF-q "^(k). By assuming A to be sufficiently large (see above), nr(k) in the first integral can be replaced by 
unity (zero) and in the second integral by zero (unity) for r = + (r = upon which one readily obtains 



^""^ The boundaries of the integral gr{<l) are determined by the observation that k in nr(k) is restricted to the range — A/2+rfcF < 
k < X/2 + rkF. 

By imposing a 'sharp' cut-off on nr(k), the long-distance behaviour of Gr coincides with that specific to FLs (Farid 1999b); 
such a behaviour is incorrect as long as a 7^ 0. 
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5r(q) = rq. 



(55) 



From this and Eq. 



it foUows that 



Urriq) 



rq 



27r iqo 



(56) 



where 11^ ^(g) stands for the density-density correlation function pertaining to the non- interacting system. This result 
is a direct manifestation of the CLT (DL 1974, Bohr 1981) which applies to the one-dimensional Luttinger model. 
Note that validity of Eq. ( [5^ ) does not depend on any other aspect of nr(k) than it be zero (unity) for r = -I- (r = — ) 
when k is sufficiently far to the right of rkp and unity (zero) for r = + {r = —) when k is sufficiently far to the left 
of rkp. Thus Eq. ( p5| ) holds also for the case where nr(k) in Eq. ( p4[ ) is replaced by its non-interacting counter-part 
0(fci? — rk); using this step function, the validity of this statement is trivially verified by explicit calculation. 

Let us now assume |q| to be small. Using (inappropriately) nr(k— q/2) — nr(k+q/2) = — q9nr(k)/5k-|-C'(|qp) for all 
k (c.f. Eq. ( p6| ) above and Eq. ( p3[ ) below), through integration-by-parts and making use of nr{X/2 + rkp — q/2) = (1) 
and n.r(— A/2 -I- rfci? -I- q/2) = 1 (0) for r = + {r = — ), up to C'(|q|'^) we evidently obtain the (exact) result in Eq. (|5q), 
in spite of the fact that nr(k) is not a differentiable function of k in a neighbourhood of k = rkp. 

We shall now demonstrate that 0(|q|^) in gr(q) — rq_ + Odq]^) is incorrect. The remainder in fact is of the form 
0(|q|^+"). Further, as the non-analytic form of this remainder suggests, the second- and higher-order expansions of 
nr(k — q/2) — nr(k -f- q/2) (in powers of q) result in divergent integrals, reflecting the fact that 1 < 1 + a < 2. 



Below we consider JX%lrk~+ii% [".(k - q/2) - n,(k + q/2)] ^ J_Z%+i/2 ^k [n.(k + rkp - q/2) - n,(k 
<^kp + q/2)] in terms of the following three constituent parts 



/•-q/2 

li:= / dk [nr(k + rkp - q/2) - nr(k -I- rkp + q/2)] , 

J-Ao/2+q/2 

/•q/2 

l2-= I dk [n^(k + rkp - q/2) - n,.(k + rkp + q/2)] , 

/■Ao/2-q/2 

Ta:^ / dk [n.r(k -I- rkp - q/2) - n^(k + rkp + q/2)] , 

Jq/2 



(57) 
(58) 
(59) 



where Ao is a constant which we chose to satisfy |q| < Ao ^ A. Our choice of the inequality Aq ^ A is motivated by the 
fact that it justifies our use of the following two-term asymptotic expression for nr(k) pertaining to the one-dimensional 
Luttinger model for spin-less fermions (Voit 1994) 



nr(k + rkp) 



■Csgn(k)|kr, 



where C is a positive constant. Making use of the expression in Eq. 
one obtains (below we assume q > 0): 



(60) 



|6C|), through exact evaluation of the integrals 



X2 



^ (Ao/2)"+i - (Ao/2 - q)"+i - q"+i 



rC 



((Ao/ 



2 - q/2)" - (q/2)" j q, 

^, (No expansion) 

(Expansion) 



2rC , 

2^-°(a+2)rC q + 1 
Q + 1 ^ 



(No expansion) 
(Expansion) 



(61) 
(62) 



The dissimilarity of the results based upon the first-order expansion of nr(k + rfci? — q/2) — nr(k+rfci? + q/2), indicated 
by "Expansion" , from the exact results, indicated by "No expansion" , confirms our above statement concerning the 
invalidity of the expansion procedure. Note for instance that T2 as determined through expansion coincides with the 
exact result only when a — 2. It is interesting also to note that contrary to the expression nj.(k — q/2) — n.r(k + q/2) = 
— q9nr(k)/9k-|-C'(|q|^) which implies absence of a zeroth-order term in q, from Eq. (60) one in fact has (c.f. Eq. (|l|) 
above) 



n^(k-q/2) - n^(k-Hq/2) 



2rCsgn(q) |k — rfc^|" -|- arC|k — vKp 

when k e [rkp — q/2, rkp 



-q/2], 



(63) 
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which contains a zeroth-order term in q as long as k ^ rkp. Without taking into account the zeroth-order term in 
Eq. (^), we would have X2 ~ 2^~"rCq""'"^ (Expansion) rather than the '(Expansion)' result for Z2 in Eq. (|62|); the 
two results coincide only for a = 1. 

Adding the three contributions Zj, j — 1,2,3, as presented in Eqs. ( |6l| ) and (|6^, we obtain 

r.AQ/2+rfeF-q/2 
'-Ao/2+rfcF+q/2 



dk [n^(k-q/2) - n^(k + q/2)] 

-An/2+rfcs-+q/2 

1^ ((Ao/2)"+i - (Ao/2 - q)"+i) , (No expansion) 



^^^i^±^q"+i + 2rC((Ao/2-q/2)"-(q/2)"). (Expansion) ^^^^ 



Assuming < q ^ Ao, and neglecting q"^^, both of the results in Eq. (|6^ ) reduce to 2^""rCAoq. This analysis 
shows that the error in (7r(q) as obtained through the first-order expansion of the integrand of (7r(q) is of order Iq]""*"^ 
and not of order |qp. As we have indicated above, in practice a is small (recall that for the one-dimensional Hubbard 
model, < a < 1/8), so that the range of the q values for which |q| is dominant with respect to |q|"''"^ is substantially 
reduced as compared with the range for which |q| is dominant with respect to |qp. Note in passing that |q|""'""'^ = |qp 
provided a — 2 which coincides with the condition required for the equality of the Z2's as evaluated without relying 
upon expansion and through (correct) expansion of the integrand of X2- 

In d > 1, contrary to c? = 1 (in particular when dealing with the one-dimensional Luttinger model), the radial 
cut-off parameter A cannot be arbitrarily large. Consequently, the condition |q| <C Aq is extremely restrictive, owing 
to Ao ^ A ^ l|kF||, and therefore should be relaxed. Correspondingly, evaluation of integrals subsequent to an 
expansion of the pertinent integrands yields results which even in the leading order deviate from the exact results (see 
Eqs. dH), (H) and (||) above). 

Finally, the momentum distribution function in Eq. (^0|) is peculiar in that nr(k) does not allow for a Taylor 
expansion in a neighbourhood of k = rfc^?. In the case of the conventional FLs, although n(k) is discontinuous at 
k = kp, nevertheless it possesses regular (Taylor) expansions to the left and to the right of k = k^?. ^ Consequently, 
in cases where one has to do with conventional FLs, provided one takes appropriate care of the discontinuity at k = kp 
of the associated n(k), the procedure based upon expansion of n(k — q/2) — n(k + q/2) yields correct results to the 
expected order. One observes that mathematical procedures that are in general inapplicable, may be applied, with 
due care, when dealing with conventional FLs (see Farid 1999a). The same procedure fails however, in the case of 
unconventional FLs or NFLs. 



V. SUMMARY AND CONCLUDING REMARKS 



In this work we have presented a brief outline of the steps that underlie a generalised method which has been 
modelled on a formalism due to DL (1974) which exactly solves the one-dimensional Luttinger (1963) model. Key 
elements in both the original DL (1974) method and its generalisation (MCDC 1998) are the conservation laws and 
the associated Ward identities. These identities give rise to some auxiliary relationships which are exact in the 
one-dimensional Luttinger model, rendering this model exactly solvable, whereas though not exact, these auxiliary 
relationships have been asserted to be asymptotically exact for Fermi systems in d > 1 in the limit where the 
scattering of the electrons in a thin layer circumscribing the Fermi surface is strictly in the forward-scattering channel 
(corresponding to scattering processes with vanishingly-small momentum transfers, that is q ^ 0). We have presented 
an analysis concerning a crucial relationship between the irreducible density and current vertex parts, A^(p;g) and 
Aa-{p; q) respectively, which is exactly valid in the case of the one-dimensional Luttinger model and asserted (MCDC 
1998) to be asymptotically exact in the forward-scattering limit when d> 1. We have explicitly demonstrated failure 
of this relationship for the conventional FL systems in d = 2. We have further demonstrated that this relationship 
is exact, in any arbitrary spatial dimension, only in the trivial case of systems of strictly non-interacting fermions. 



Within the framework of the RPA, n(k) pertaining to the uniform system of electrons (in d = 3) has been shown to possess 
finite left and right derivatives at k = k_F (Daniel and Vosko 1960). On the other hand, for a uniform system of electrons 
interacting through a hard-core two-body potential, Belyakov (1961) (see also Sartor and Mahaux 1980) has obtained an n(k) 
which has logarithmically divergent derivatives to the left and right of k = k^. For a discussion of the behaviour of n(k) in the 
close vicinity of k = k_F see (Farid 1999a). 
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From the same analysis it follows that for NFLs in d > 1 there does not exist even the remotest ground on the 
basis of which the presumed asymptotic validity of the mentioned relationship can be justified: here, the interlocking 
of energies and momenta in the arguments of the corresponding self-energies, which by the very fact of being by 
definition non-differentiable in neighbourhoods of the Fermi energy and/or the Fermi momentum (otherwise the 
associated metallic states would not be NFLs) , stands in the way of obtaining a relationship similar in appearance to 
the indicated supposed asymptotic relationship between the irreducible vertex parts. 

We have further explicitly demonstrated that the continuity equation only in d = 1 gives rise to a direct relationship 
between the vertex parts AJ^(p; q) and r[^(p; g), i/ = 0, 1, . . . , d, independent of the analytic properties of the momentum 
distribution function nCT(k). In d > 1, continuity as well as continuous differentiability of this function in the above- 
indicated narrow band of the momentum space circumscribing the Fermi surface must be assumed before a similar 
direct relationship between the two vertex parts, as in d = 1, can be deduced. As is well known, no-(k) is discontinuous 
at k = both for FLs and certain class of NFLs (for the latter see Farid 1999a). For the case of the one-dimensional 
Luttinger model, though continuous, nCT(k) is not differentiable at k = \ip. Therefore, by assuming a similar behaviour 
for ncr(k) in d > 1 as in d = 1, we have established that the indicated simplified relationships between AJ^(p;(j) and 
T'^{p; q), 1/ = 0, 1, . . . , d, do not hold in d > 1. 

Each of the above-mentioned two findings removes ground for validity in d > 1 of the heretofore supposed asymptotic 
Ward identities in the forward-scattering limit. This failure concerns both FLs (excluding strictly non- interacting 
metallic systems) and NFLs. 

In an earlier work (Farid 1999b) we have demonstrated failure of the Haldane (1992) bosonisation scheme in d > 1 
both on general grounds and on grounds of an explicit calculation on a model of interacting fermions in d = 2. Within 
the framework of the Haldane bosonisation scheme, the explicit expression for the single-particle Green function G 
takes a simple form in the space-time representation. It is interesting that through replacing a screening function that 
features in this expression by its random-phase approximation (that is RPA) counterpart, one obtains an expression 
which coincides with that obtained within the generalised DL framework; compare Eq. (5.65) in (MCDC 1998) with 
for instance Eq. (24) in (Kwon, Houghton and Marston 1995). The use of the RPA within the framework of the 
generalised DL method is considered to be justified through reliance upon the so-called generalised CLT (Kopietz, 
Hermisson and Schonhammer 1996) which implies that while the RPA yields the exact response function in the case 
of the one-dimensional Luttinger model (DL 1974) (see Eq. ( [s^ ) above), it yields an asymptotically exact response 
function in d > 1 in the forward-scattering limit. In (Farid 1999b) we have commented on the generalised CLT, 
however independent of the true status of this theorem in d > 1, and whether one has to do with FL or NFL systems, 
from the point of view of our considerations in the present work the coincidence of the two Green functions is of crucial 
importance as it indirectly establishes that indeed the generalised DL (MCDC 1998) formalism is of no validity for 
systems of interacting fermions in d > 1. 

To summarise the main conclusions arrived at in this and our earlier works (Farid 1999a, b), we can state that 
none of the presumed non-perturbative methods designed thus far for investigating the existence or otherwise of NFL 
metallic states in d > 1 has proved capable of dealing with systems of interacting fermion in d > 1, not even in 
the weakly- interacting limit. Renewed attempts concerning design and application of non-perturbative methods are 
therefore called for in order to settle the long-standing problem with regard to the existence or otherwise of NFL 
metallic states in spatial dimensions greater than one. It is therefore befitting to close this work by mentioning that 
Laughlin's (1998) recent statement, that "... and all attempts to account for the existence of a Luttinger-liquid at zero 
temperature in spatial dimension greater than one have failed." , must needs be qualified and in due course perhaps 
revoked altogether. 

□ 
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Recall that metallic states of non-interacting fermions are ideal FLs. 
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APPENDIX A: ON DIMENSIONAL CONTINUATION AND ITS APPLICATION TO INTERACTING 

FERMI SYSTEMS 

In this Appendix we first present a brief account of the technique of dimensional continuation (Wilson 1973, 
appendix). ^ Subsequently, we point out that the formalism of interacting fcrmions even in the limit of strong 
forward scattering cannot be subjected to the conventional process of dimensional continuation. We consequently 
propose an unconventional approach to this process which involves one continuous free parameter and therefore 
defines an infinite family of continuations that correctly interpolate between d — 1 and d — 2; they also correctly 
interpolate between d — 2 and d = 3 for the class of cylindrically symmetric integrands. 

Consider 



Jrf(q):- / d'^k/q(k). 



(Al) 



where d is some real or complex number. Let {ei, . . . ,ej} be an orthonormal basis spanning the so-called parallel 
space (to whose complimentary space we refer as tangent space), large enough to accommodate the vector q (in 
general, all external momenta {qj} upon which the integrand may depend). We have 



S = 5|| U5j.. 



(A2) 



With 



k = k|| + kx = ^ fc|ej + k^, 
i=i 

one can therefore write 

Jd(q)= / d^...dfc^ I d'^-^ki /q(k). 

5 II >/ S'_i_ 

By assuming independence of /q(k) upon direction of k^, i.e. by assuming 



(A3) 



(A4) 



/q(k) =/q(k.k,k-q) 



(A5) 



one has 



Js, Jo 



dllk^ll l|kj 



Id- J- 



i/q(kk,k.q), 



(A6) 



where 



(A7) 



with r(z) the gamma function. The value of K^, in Eq. (A_7) is equal to the area of the hyper-sphere of unit radius 
in I' dimensions. The J + 1 integrals over fc||, . . . ,k^l and ||kj^|| are conventional and are evaluated in the usual way. 



It should be evident that the key element for reducing the abstract integral in Eq. (Al) into a conventional multi- 
dimensional integral is the assumption specified by the expression in Eq. (A5). Below we explicitly consider two cases 
which clarify the role played by the dimension of the parallel space, J. 



Case I: Let J = 1. Let further 6 denote the planar angle between k and q. With 



For a comprehensive review of this topic see Collins (1984, Chapter 4). 
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||k_L|| = I|kl|sin6l, fc^ = ||k||cos6', 6lG[0,7r], (A8) 
one has dfc^d||ki|| = ded\\k\\ Hkj^f-^ ^ ||kf sin6')''"^ and thus (Wilson 1973) 

/OO nOO 
dklKd-i / dijkxil |!kxf-'/q(k.k,k.q) 
-OO Jo 

= Kd-i d9 d||k|| ||kf-i(sin0)''"'/q(||k||M|k||||q||cos0). (A9) 
Jo Jo 



Since Kd^i{smd)'^ ^ 5{d) + S{9 — vr) for d J, 1, it follows that in the limit d i 1 the expression on the RHS of 



Eq. (A£) recovers the pertinent expression for d = I, i.e. 

JM)\dll = / d||k|| /qdlkf , ||k|| llqll) + / d||k|| /q(||kf , -||k|| ||q||) 
Jo JO 

dfc/q(fc2,A:||q||). (AlO) 



For d — 2, for which Kd-i = 2 holds, we have 

nTT nOQ 

M^)\d=2='2 dW d|[k|M|k||/q(||kf,||k||||q||cos0) 
JO Jo 

^ d^ d||k|| ||k||/q(||k||M|k||||q||cos^). (All) 
JO JO 

The latter result, whose validity can be readily verified, |^ makes explicit that for d = 2 one indeed recovers the 
conventional integral over k of /q(||k|p, k • q) in terms of the cylindrical-polar coordinates (||k||, (^), with q the polar 
axis. It should be noted however, that not all two-dimensional integrals can be viewed as the d 1 2 limit of the 



expression on the RHS of Eq. ( A9) 



Case II: Let J = 2, so that {61,62} spans the parallel space. With Lp denoting the planar angle between k|| and ei 
(which without loss of generality we assume to be along q) and 6 that between k and kj^ (i.e. kj_ is taken as the 
polar axis), one has 

|lk_L|l = ||k||cos6', fc^ = ||k||sin6l cos^, fcjj = ||k|| sin6l sin(^, 6'G[0,7r/2], <^e[0,27r], (A12) 
so that dfc^dfcjjd||ki|| = d^jdMllkll ||k||2sin6', Hk^f-^ = ||kf ^^g^-j'^-s ^^^^ 



/OO noc 
dk\dklKd-2 / dllk^ll ||kif-Vq(k-k,k.q) 
-00 Jo 

/ d^ / de'(cose') sin6l 



Jo 

X 



poo 

/ d||k|| ||kr-i/q(||k||M|k||||q||sin0cos^). (A13) 
Jo 

Since Kd-2{coaey^^ -> 6{e - tt/2) for d i 2, it follows that 

p27T poo 

JM)\di2^ d^ / d||k||||k||/q(||k||M|k||||q||cos^), (A14) 
Jo Jo 



To this end employ the variable transformation if ^ u defined through u = cos 75, implying d;^ = =pdu/(l — u^Y^^ for 
if G [0, tt] and if G [vr, 27r] respectively, and a subsequent reverse transformation u ^ f. 
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which is indeed the expected expression for a two-dimensional integral in the cylindrical-polar coordinates (||k||,(^). 
However, again as in the case considered above (concerning J = 1), not all two-dimensional integrals over arbitrary 
functions of k can be obtained as the d i 2 limit of the d-dimensional integral in Eq. ( A13). 



The above considerations make explicit that whereas by an appropriate choice of the dimension of the parallel space 
one can analytically continue d and obtain the standard integrals for integer dimensions, nonetheless the applicability 
of this procedure in an essential way depends on the requirement that /q(k) in so far as its dependence upon k is 
concerned, be a function of at most k • k and k • q. In cases where the dependence of /q(k) upon k is more general, 
the conventional process of dimensional continuation may be possible in the limit of small l|ql|. In this limit, provided 
/q(k) be continuously differentiable in a neighbourhood of q = 0, it may be possible (that is, under the condition that 
/q(k) to linear order in q can be described by a function of k • k and k • q) to carry through the process of dimensional 
continuation, accurate to linear order in q. ^ This is precisely the case in the generalised DL formalism by MCDC 
(1998) (see the last part of § II). 



It is readily verified that the following expression, whose origin does not lie in the conventional formalism for 
dimensional continuation (Wilson 1973), is a smooth analytic continuation for d G (1, 2]: 



Jd(q):= 



K, 



d-i 



m wo 



d0 (sin^ 



, d-2 



d9 (sin(6l-7r)) 





d-2 



d||k|| ||kr-l/q(||k|| 



'd||k|| ||kr-l/q(||k||,0)}. 



(A15) 



where 



/q(||k||,0)^/q(k), 

with (||k||,6') the cylindrical-polar coordinates of k, and (see Figure) 

^9(d):=l + (3-d)exp(a[l-l/((i-l)i/'']), with a > 0. 



(A16) 



(A17) 



We note that the integral in Eq. ( A15| ) together with 'd{d) in Eq. (A17) is also correct for d = 3 over the space of 
functions with cylindrical symmetry; therefore it is also a smooth continuation for d £ [2, 3], however over a restricted 
set of functions. In this respect Eq. ( Al5| ) is of limited validity for d > 2, in the same way that Eq. (A9) is of limited 

validity for d > 1. 

The choice for ?9((i) in Eq. ( Al7| ) is not unique. However, in choosing any 'd{d) interpolating between ?9(1) = 1 
and 13(2) — 2 (and possibly 19(3) = 1, etc.) care must be exercised that in the limit d 1 the singularities of 
Kd-i{smd)'^ ^ a.nd Kd-i{sm{6 ~ tt)^"^ ^ 
is prevented by the requirement {'d{d) 



I TT and J, TT do not yield repeated contributions. This 



1 



for respectively 

^'^ ^ ^ for d I 1; in consequence of this, for d I 1 the reduction in the 

width of the interval [TT,d{d)Tr] renders the growth in the integrand, due to the singularity of Kd-i(^sm{9 — tt))'^ ^ 
at d = 1, ineffective. It can be easily verified that 'd{d) in Eq. ( |A17 ) indeed satisfies this requirement. The freedom 
in choice of the continuous positive parameter a in Eq. (A17) implies that the dimensional continuation scheme as 
prescribed by Eqs. (A15) and (|A17 ) represents an infinite family of dimensional continuations. However, in order to 
recover the expression in Eq. ( |A9D , applicable to functions {/q(k)} whose dependence upon k is through k • k and 
k • q, it is required that < a <C 1 (see Figure); in the limit a I the expressions in Eqs. (|A9|) and (A15) coincide. □ 
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